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We investigate the efficiency of cooling the vibrations of a nano-mechanical resonator, constituted 
by a partially suspended Carbon- nanotube and operating as double-quantum dot. The motion 
is brought to lower temperatures by tailoring the energy exchange via electromechanical coupling 
with single electrons, constantly flowing through the nanotube when a constant potential difference 
is applied at its extremes in the Coulomb-blockade regime. Ground-state cooling is possible at 
sufficiently high quality factors, provided that the dephasing rate of electron transport within the 
double dot does not exceed the resonator frequency. For large values of the dephasing rates cooling 
can still be achieved by appropriately setting the tunable parameters. 



I. INTRODUCTION 

The experimental demonstration of strong coupling 
between a nano-mechan ical resonator and charge trans- 
port in carbon nanotubes^' holds promising perspectives 
for cooling the temperature of the resonator vibrations 
by means of electro-mechanical forces. It may enable 
one to improve the sensitivity of high-precision measure- 
ments of mass^H^, mechanical displacement, or spin ^. 
Cooling by means of electro-mechanical coupling, more- 
over, may allow one to prepare the resonator state in 
the quantum regime 3 -'. Recently, it has been theoreti- 
cally proposecP^HUl and experimentally demonstrated^ 
that back-action cooling can be achieved by coupling 
mechanical resonators to the constant electron current 
flowing through electronic nano-devices, such as normal- 
metal and superconducting single-electron transistors^. 
Within these studies, it is found that one of the major 
factors limiting the cooling efficiency is the quality factor 
of the mechanical resonator. In this respect a suspended 
carbon nanotube is a promising candidate for achieving 
ground state coohng-221. 

In this work we theoretically analyze the efficiency of 
ground-state cooling the vibrations of a nano-mechanical 
resonator using constant electron current, in a setup like 
the one sketched in Fig. [T] Here, the resonator is a sus- 
pended nanotube, which constitutes also the electronic 
device through which the current flows, in a s etup which 
is reminescent to the ones realized in RefsP^l. Moreover, 
the nanotube acts as a double quantum dot (DQD). This 
setup allows us to tune the electron current, so to enhance 
the electromechanical processes which irreversibly absorb 
phonons from the resonator. The analysis here presented 
extends and complements the proposal in Ref!^, where 
we predicted ground-state cooling in this system for a 
specific set of parameters. In the present work we study 
the cooling efficiency over a larger parameter space than 
the one identified in RefM^U, showing that large ground- 
state occupations can also be obtained for relatively large 
tunneling rates. Moreover, we analyze in detail the role 




FIG. 1: (Color online) (a) Sketch of the considered set-up. 
The double quantum dot is a carbon nanotube (CNT) op- 
erating as single-electron transistor. In the picture, the dot 
on the right is suspended and mechanically vibrates. Its vi- 
brations couple to the electronic current via electron-phonon 
interaction. In this article we will also consider the situation 
in which it is the dot on the left which is suspended, (b) 
Energy diagram for the DQD. States \L) and \R) denote one 
excess electron in the left and right dot, respectively. The 
two levels are at tunable energy difference hA, and are cou- 
pled by a tunneling barrier with tunneling matrix element T c . 
Left (right) electrode act as source (drain) such that electron 
tunnels from left to right at rates Tl and F_r. 



of the various noise sources on the cooling efficiency. 

It must be mentioned that cooling and manipulation 
of nano-mechanical resonators using electro-mechanical 
forces is complementary to the approach using the cou- 
pling to photoniP^Hall Both approaches are presently ex- 
perimentally pursued. In particular, cooling by means of 
constant electron currents is appealing because it is easy 
to implement in a dilution fridge as compared to tech- 
niques based on photons. Beyond the specific implemen- 
tation, ground-state cooling of the mechanical resonator 
would open the possibility to create and manipulate non 
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classical states at the macroscopic scale and to st udy the 
transition from the classical to the quantum regim o 3 l 32 l 33 [ 
This article is organized as follows. In Sec.|Tl]the model 
is introduced, which describes the coherent and incoher- 
ent dynamics of the nanotube. The rate equations for 
cooling are derived in Sec. and the cooling efficiency 
is discussed as a function of the physical parameters. The 
conclusions are presented in Sec. |IV| 



with 



II. THE MODEL 

A. Setup and basic dynamics 

The resonator we consider corresponds to a fundamen- 
tal vibrational eigenmode of a partially suspended carbon 
nanotube (CNT). The eigenfrequency of the mode is lj, 
and the resonator is assumed to be characterized by suffi- 
ciently high quality factors, so that the eigenmode can be 
selectively addressed by suitably tailoring the mechani- 
cal coupling with the constant electronic current flowing 
though the CNT. The eigenmode is here described by a 
quantum harmonic oscillator, with a' and a the bosonic 
creation and annihilation operators of a quantum of en- 
ergy Tiuj. The Hamiltonian for the resonator is given by 

H p h — fiuja) a , 

and it governs its coherent dynamics. The eigenstates 
are the number states \n), with energy Timo and n = 
0,1,2,... 

Before discussing the details of the electromechanical 
coupling, we first introduce the electronic properties of 
the system. The system is in a double quantum dot 
(DQD) configuration, as sketched in Fig ija): Voltage 
applied to gate T provides the potential barrier, which 
divides the CNT in two coupled quantum dotd 3 -^ 3 -^. The 
voltage applied at the left and right electrodes warrants 
that a constant electron current flows through the CNT 
from left to right. We denote by /i£ and n R the chemi- 
cal potential at the left and right electrode, respectively, 
such that hl > fi R . We assume that the DQD is working 
in the Coulomb-blockade regime such that the charging 
energy E c is much larger then the transport energy win- 
dow, E c ^> — fiR. In this regime the system operates 
as a single-electron transistor (SET), and the relevant 
electronic states are the ground state of the DQD, |0), 
and one excess electron either in the left or in the right 
dot, corresponding to the quantum states \L) and \R) , re- 
spectively. More specifically, the electronic ground state 
corresponds to the state with Nl and N R electrons in 
the left and right QD respectively, |0) = \Nl,Nr), while 
\L) = \N L +_LNr) and \R) = \N L ,N R + 1), see for in- 
stance RefspSHini We denote by el and e R the energy of 
the states \L) and \R) with respect to the ground state. 
Their difference is controlled by the gates L and R in 
Pig Qa). The DQD Hamiltonian, governing the dynam- 
ics of the electrons in the DQD, reads 



H LR = e L \L)(L\+e R \R)(R\, 
H T = -hT c (\R)(L\ + \L)(R\) 



(2) 
(3) 



where T c is the tunneling matrix element between left 
and right QD and has the dimension of a frequency. The 
energy levels of the QDs and the direction of the electron 
current are depicted in Fig[l|b). 

We now discuss the electromechanical coupling. The 
vibration are capacitively coupled to the gate R: The gate 
voltage bends the CNT influencing its electric and me- 
chanical properties. When the vibration amplitude is suf- 
ficiently small, the electromechanical energy can be con- 
sidered to be linear in the amplitude. Assuming that the 
energy difference between left and right dot is such that 
e R — €l ~ Twj, the electronic transport couples with the 
fundamental eigenmode at frequency w, and this mode 
is displaced along x so that the Hamiltonian describing 
this dynamics can be approximated by the operator 



H e _ ph = ha\R){R\ (a^+a) 



(4) 



where we used x oc a) + a, and a has the dimension of 
a frequency and gives the strength of the electromechan- 
ical coupling. This Hamiltonian describes a shift of the 
center of oscillation of the harmonic vibrations, condi- 
tioned to the occupation of the electronic state \R) of 
the DQD or, in other terms, an energy shift of the state 
\R) proportional to the mechanical displacement of the 
resonator. 

We finally introduce the model for the coupling of the 
DQD with the electrodes, which gives rise to the constant 
electron current flowing through the CNT. The electrodes 
are electron reservoirs, described by the Hamiltonian 



k k 



(5) 



where and dk (cl and dV) are fermionic operators, 
which annihilate (create) an electron of the respective 
reservoir and k runs over all modes of the electrons. Elec- 
trodes and QDs are coupled by electron tunneling, which 
in Hamiltonian form reads 



W L = E^l^ + H.c. 

k 

w r = Y. v * d i s L+K. c . 



(6) 
(7) 



H. 



DQD 



Hlr + Hj 



(1) 



where V k L and V k R are the tunneling matrix elements be- 
tween reservoirs and quantum dots, and = |0)(j| an- 
nihilates an electron in the QD j, with j = L, R, leaving 
the QD in the ground state while an electron is created 
in one mode of the reservoir. 

In what follows we assume that the bias voltage is suf- 
ficiently large such that thermal effects in the coupling 
with the electrodes are negligible. This implies that the 
relations /i^ — fi R 3> ksT and [1l 3> £l, £r 3> must 
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hold. In this regime the electrons tunnel from left to 
right, while transport in the opposite direction can be 
neglected. Tunneling from and into the electrodes can be 
described in perturbation theory. In first-order perturba- 
tion theory the tunneling rate T l from the left electrode 
into the DQD, and the tunneling rate Tr from the DQD 
to the right electrode, are given by 



R\ 
k J ! 



(8) 
(9) 



where <5(e) expresses the condition of energy conservation 
between initial and final state. In the continuum limit the 
sums over the modes of the electrodes become integrals, 
properly weighted by the density of states. 



B. Master equation 

In order to describe the effect of quantum noise on the 
resonator dynamics we resort to the density matrix for 
the degrees of freedom of the mechanical resonator and 
of the DQD electronic degrees of freedom. Its dynamics 
is governed by the master equation, which also accounts 
for incoherent tunneling from and into the electrodes, 
dephasing of trasport between the QDs, and thermaliza- 
tion of the resonator to the temperature at the electrodes. 
The master equation reads 

d 1 

-Q t P = [Hdqd + H ph + H e - ph , p] (10) 
+£dqdP + Cdp + ICp , 

where the term 

£dqdP = T L /2 (2s\ps L - ps L s\ - slS^p^J 

+T R /2 (2s R ps R - ps R s R - s R s R pj (11) 

describes incoherent electron tunneling from left elec- 
trode to state |L), and from state \R) to the right elec- 
trode, while the term 



C d = l -f(2 S ^p S ^-s^ 2 p-p S ^ 2 ) 



(12) 



describes electronic dephasing inside the DQD at rate Yd, 
with sW = \R)(R\ - \L)(L\. Finally 

JCp = (hp + l)7 p /2 [2apa i — a^ap — pa) a) 

+h p "fp/2 (2a' pa — aa' p — paa') (13) 

represents processes leading to thcrmalization at rate j p 
with a thermal environment at temperature T, such that 
the mean number of phonon at frequency lo is given by 

hp = [exp (fioj/kBT) — 1] 1 . 



The Liouvillian terms in Eq. ( 10 1 are written in the Born- 
Markov approximation: The master equation is local in 
time, it hence does not consider memory effects which 
may arise from the coupling with the reservoir. A mas- 
ter equation for a double-quantu m dot coupled to a res- 
onator has been reported in RefsPsEH. With respect to 
that master equation, here we have added the finite tem- 
perature of the electrodes and the dephasing mechanism. 
The latter is here modeled according to the master equa- 
tion presented for instance in RefP^l. 



III. ELECTRO-MECHANICAL COOLING 

Ground-state cooling via electro-mechanical coupling 
can be achieved by enhancing the processes of electron 
transport which absorb a phonon of the nano-mechanical 
resonator. This is realized by properly setting the energy 
difference between the two QDs. Such processes must 
take place at rates which overcome the rate of electro- 
mechanical processes which tend to heat the resonator 
and the rate at which the resonator naturally thermal- 
izes with the surrounding environment. In this section 
we derive the basic equations which describe the cooling 
dynamics of the resonator, and study them in different 
limits. The validity of the analytical predictions is tested 
by comparing them with the numerical solution of the 
master equation. 



A. Rate Equations for cooling dynamics 

In order to identify the important parameters, which 
determine the cooling dynamics, we resort to perturba- 
tion theory and derive rate equations for the occupation 
of the resonator's vibrational states. This is performed 
eliminating the electronic degrees of freedom from the 
resonator dynamics, assuming that a, scaling the elec- 
tromechanical coupling, is much smaller than the char- 
acteristic rates determining the electronic dynamics. In 
this regime, we assume that the electronic degrees of free- 
dom of the DQD are in the steady state. This steady 
state determines the properties of the electronic cur- 
rent, and thus also the autocorrelation spectrum of the 
electro-mechanical force acting on the resonator. Assum- 
ing time-scale separation between the electronic dynam- 
ics and the dynamics of electromechanical coupling, we 
then derive a master equation for the resonator degrees 
of freedom, obtained from master equation ( 10 1 in second 



order perturbation theory and tracing out the electronic 
variables. For this purpose we rewrite the master equa- 



tion ( 10 ) as 



0_ 

Wt 



p = (C a + Ci + K.) p, 



where 



C p = Co 'p 



[H p h,p] , 



(14) 



(15) 
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accounts for DQD and resonator dynamics in absence of 
electro-mechanical coupling, with 

£ ° ) p = -^[H DQD ,p] + (£ mJ3 +£ d )p (16) 

the Liouvillian for the DQD. The term £\ contains the 
electromechanical coupling term, 



£\P 



-[H e - P h,p] , 



(17) 



and it scales with a. Finally, the Liouvillian K, describes 
the thermalization of the resonator, it scales with 7 P , and 
is given in Eq. (13 1. 



We assume that the back-action of the electromechan- 
ical coupling on the electronic state can be neglected, 
which corresponds to considering the regime where |a| <C 
Tr,Tl,\T c \. Then, the electronic steady state pg°' is de- 
termined by the solution of the equation 



£ 



(o) Jo) _ 
o Pst — 



0. 



(18) 



In particular, p^ can be written as 

p$ = p 00 \0)(0\+p LL \L)(L\+p RR \R)(R\ 
+PB. L \R)(L\+p LR \L)(R\ 



(19) 



where the elements scaling with the parameters pg R and 
Pol, and corresponding to the coherences between vac- 
uum and \L), \R), are not reported, as /?ol and po R vanish 
at steady state due to incoherent tunneling between the 
electrodes and the DQD. The other parameters take the 
form 



1 



Pll = 



£2 

4 



Prr 



Prl 



= 2- 



A 2 



TV rp2 



2 £^ + r' 2 



£2 

4 



(A + ir^) 



(20) 
(21) 
(22) 



with p LR = p* RL and p 00 = 1 - p RR 



(22) we defined T = 2 



r' 



lit 



2T; 
T7 



Pll- In Eqs. (20) 

1/2 



r' 



and 



r fl /2, 



(23) 



giving the total decay of the electronic coherences of the 
DQD. Wc further assume that 

\a\ <w. 

In this limit, a closed set of equations for the occupation 
of the resonator states can be derived. For this purpose, 
we define the projector P acting over the density matrix 
of the system, such that its action over a given density 
matrix p of the system reads 



Pp 



(o) 
Pst 



Tr cl {P p} • 



(24) 



where Tr c i denotes the trace over the electronic degrees of 
freedom and Pop — J2 n \ n )( n \( n \p\ n )- The master equa- 
tion for the density operat or pro jected over this subspace 
can be written in the for m 1 43 * 44 -! 



Pp = PJCPp + P / dr£ ie £ ° T £ 1 Pp(t), (25) 



d_ 
Wt 

where we applied a perturbation expansion to second or- 
der in a, and took 7 P ~ o(a 2 ), so that thermalization 
effects arc taken at lowest order in the expansion. This 
latter assumption is indeed important for the efficiency of 
the cooling process, as only in this limit the electrome- 
chanical coupling may be able to counteract effectively 
thermalization with the external environment, bringing 
the resonator into a dynamical equilibrium with the elec- 
tronic current. 

After tracing out the electronic degrees of freedom 
from Eq. ( 25 ) one obtains a set of rate equations for the 
occupation p n of the reson ator state with n phononic ex- 
citations, which rea d 45 * 46 ! 

Pn = (n + l)A_p n+1 - [(n + l)A + + nA-} p n + nA|_p„_i , 



(26) 



where A + and A- are the rate for processes which in- 
crease and decrease, respectively, the state of the me- 
chanical resonator by one excitation. They result from 
the interplay between the thermalization processes and 
of the cooling due to the mechanical effects induced by 
the electron current. Within the parameter regimes con- 
sidered so far, they are the sum of thermalization and 
electro-mechanical rates, and written as 



-44 



1) 



A -, 



(27) 

i+ = 7 P n p + A 0+ , (28) 

where A Q+ = 2Re{S{-u)} (A _ = 2Re{S(u)}) is the 
rate for electromechanical heating (cooling), with 



S(u J ) = -a 2 Tr\\R)(R\ 



the spectrum of the autocor relation function of the elec- 
tromechanical coupling forc d 43 * 44 ^. It can be rewritten 

as 



5(w) = p RR S RR {u) + prlS rl (u) , 



(30) 



where 
Srr(uj) 

Srl(oj) 



2F(uj) 
. a 2 T, 



(r ^- iw ) 2 + 2i^(r^-i w )4-A 2 



2F(u) 



[r' d -i( w + A)] , 



(31) 



and F(ui) = ^fF R - i^Fi, with 



Fn 



r 2 

T 



cu 2 I 1 



Ft = 2t;t 



A 2 

Tl - n 



2r^ , 2T 2 r^-r, 

Ifl 



Tl r 2 



(T / d + 2T R ) 



+4T 2 



(32) 
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In the following analysis we will also consider the situ- 
ation, in which the resonator couples only with the left 
dot. In this case, the equations are modified accordingly, 
starting from the spectrum of the autocorrelation func- 
tion of the electromechanical force, which in this latter 
case reads 

S'(u>) = -a 2 Tr{|L)(L| (4» +iu,) _1 \L)(L\($ 

The corresponding analytical form of the rates are not 
reported here, but are derived using the procedure de- 
scribed above. 



B. Cooling dynamics 

From the rate equation for the population of the 
phononic states one simply finds an equation for the 
mean num ber o f mechanical excitations n = "}2 n np n , 
which readsP^ 21 



the average excitation number at steady state, 



-7tot« + A + , 



(34) 



where 7tot = A- — A + is the cooling rate, and in partic- 
ular 



7tot = 7 P + 7o 



(35) 



with 7q = A _ — A 0+ the electromechanical cooling rate. 
A steady-state solution is found when 7 tot > 0, giving 



"-st 



n P 7 P + n 7o 
7 P + 7o 



(36) 



where no = ^4-o+/7o is the mean phononic number when 
7p is set to zero. The cooling dynamics results from the 
competition between the mechanical effect of the elec- 
tron current and the thermalization with the environ- 
ment. Clearly, the electromechanical coupling cools the 
resonator below the temperature T of the environment 
when ngt < n p , which is possible when no < n p . Equa- 



tion ( 36 1 shows that lowest temperatures are achieved 



when 70 > 7 P , giving 



"st ^ n + (rip - n ) 7 P /7o ■ 



In the following we will search for the parameters which 
minimize fist for different values of the thermalization 
rate 7 P and of the dephasing rate IV In order to maxi- 
mize the cooling rate, we first assume that Tl, the rate at 
which electrons are injected into the DQD from the left 
electrode, is the largest rate characterizing the dynamics 
(the higher bound is provided by the condition that the 
system must operate in the SET regime), implying that 
the time intervals, in which there is no excess electron 
inside the DQD, are here the smallest time scales. In the 
limit of large injection rates, the cooling rate 70 reads 



3 A 2T 2 



7o 



A 2 + 2T 2 



O 



(A 2 + 4T 2 



4 1 R 



r 2 

1 R 



A 2 + 2T 2 



2w 2 



. (37) 



while the mean phonon number n takes the form 



n 



r 2 ? + 4(A-c) 2 

16Aw 



o 



(38) 



These expressions have been obtained setting the dephas- 
ing rate Td -C Tr. The effect of dephasing rates Td ~ 
on the cooling efficiency will be discussed later on. These 
expressions are valid both for the case in which the res- 
onator couples to the right or to the left dot. Differences 
in the cooling efficiencies for the two setups arise when 
the value of Tl is lowered, so that it becomes comparable 
with other rates, as we will show. 

Ground-state cooling. Ground-state cooling may be 
obtained when fio <SC 1 and 70 3> 7 P . Assuming large 
injection rates into the DQD, and setting the dephasing 
rate Td = 0, small values no <§C 1 are achieved in Eq. (38) 



for uj 3> Tjj. The minimum value 
no = r|j/I6w 2 



is found setting A = oj. From Eq. (37) we find that 
the cooling rate is maximum when tunneling rate and 
detuning fulfill the relation lo = e with 



e = v/A 2 + 4T 2 . 

This condition corresponds to the one which minimizes 
the temperature, provided that \T C \ <C uj. At suffi- 
cienctly large tunneling rates this condition emerges from 
the physical situation, in which the electronic states of 
the DQD are described by the bonding and antibond- 
ing states, such that e is their frequency splitting. These 
states are the eigenstates of the Hamiltonian 0, and 
are given by the quantum superposition of left and right 
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electronic states 



cos6\L) +sm0\R) 
-sin9\L) + sin6»|i?) 



(39) 



with tan 8 — 2T C / (A + e). In this basis there are two 
relevant physical processes, which lead to a change by 
one phonon due to electron transport through the DQD. 
They consist of the sequential occupation of the states 
|0, n) -> \±, n) — ¥ n =F 1) -> |0, n^l). Both processes 
are resonant when the condition e = uj is satisfied. The 
resonator is cooled when the rate of the cooling process 
is faster than the heating process, which is here satisfied 
for A > 0, see Eq. ([37} and RefP. 

Figure [2] displays the contour plots for fist and 7 tot 
for a mechanical resonator with quality factor Q = 10 5 , 
the results are compared with the ideal situation j p — 
(Q —> oo). We first discuss the ideal case. Here, we 
notice that the final mean occupation depends solely on 
the detuning A, as expected from Eq. ( [38] ), while large 
cooling rates are achieved when T c and A satisfy the res- 
onance condition uj 2 ~ A 2 + 4T 2 and < A < u. At 
finite Q, efficient cooling is found for the condition at 
which the cooling rate is maximized, as expected from 
Eq. ( 36 1 . Note that the cooling efficiency decreases both 
in the limit T c — >• (A — >• to) and A -> 0. In the first 
case, the total electron transport rate decreases, and cor- 
respondingly also the rate of electromechanical processes. 
In the second case, the asymmetry between heating and 
cooling rates (A±) is reduced. 

Effect of dephasing on the cooling efficiency. We now 
proceed in studying the effect of dephasing on the ef- 
ficiency of the cooling mechanism. In Fig. [3] average 
phonon number and cooling rate are reported for Q = 10 5 
and setting Yd = O.lw, 5ui. Clearly, the cooling efficiency 
is lowered as the dephasing is increased. In particular, 
the larger is the dephasing, the slower becomes the cool- 
ing. Nevertheless, for = 5uj it is found that the res- 
onator can be cooled from n p — 50 to n = 10 with a 
rate 7 tot = 10~ 4 w. We also observe that the parameter 
region, where cooling is achieved, is reduced (recall that 
the initial phonon number here considered is n p = 50). 
Optimal cooling is obtained for values of the detuning A 
of the order of T^. Indeed, a simple calculation of the 
average phonon number, taking Yd > us,Tr,T c , shows 
that this is modified according to the formula 



A 2 



2Aw 



(40) 



which reaches the minimal value n' Q ~ T^/o; at A = T^. 
Correspondingly, the cooling rate is larger when to > Tr 
and scales with 70 ~ a 2 T 2 / (T 2 d uj) . The cooling rate hence 
increases with the ratio between coherent tunneling and 
dephasing rate, with the upper bound 70 < a 2 /uj. This 
cooling limit is reminiscent of Dopplcr cooling of atoms 
in a dipolar transition with effective lincwidth 2Td and 
recoil frequency ujr ~ a 2 /uj, see Refs.™ 

Cooling efficiency as a function of the incoherent tun- 
neling rates. So far we have assumed that injection rate 
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FIG. 2: Contour plots of the average phonon number at 
steady state fist, Eq. ( |36[ ), and of the cooling rate jtot, 
Eq. (341 as a function of A and T c , in absence of dephasing. 



The grey scala is such, that the darkest (brightest) regions in- 
dicate the smallest (largest) values, with the exception of the 
region labeled with "no steady state", where the resonator 
is heated by electromechanical processes. The dashed line in 
the plots (b) and (d) indicates the curve 

A 2 + 4T 2 _ (J 2_ The 

parameters are Tl = IOOoj, Tr — 0.2cj, a = O.Ioj, 1^ = and 
(a)-(b) 7 P = 0, (c)-(d) 7 P = 10" 5 cj with n p 



50. 



into the DQD, Tl, is the largest parameter, while we 
took < uj. We now study the cooling efficiency as a 
function of and T^. We first focus on the situation in 
which is the largest parameter and consider the cool- 
ing rate as a function of T^, when T' d ~ Tu/2. Figured 
displays the average excitation number nst and the cooh 
ing rate 7tot as a function of T^, choosing T c and A so 
to optimize both parameters. One clearly observes that 
optimal cooling is obtained at small values of T^, while 
the cooling rate vanishes at Tr — > 0. The numerical sim- 
ulations show that the cooling efficiency has a maximum 
at Tr of the order of a fraction of the oscillator frequency 
uj, while it decreases as Tr is further increased. We note, 
however, that phononic occupations fist < 1 are still ob- 
tained for ~ uj. The analytical results reproduce well 
the behaviour at larger values of r#, where one finds 
70 ~ r^j 2 . A significant discrepancy between analytical 
and numerical results is observed in the limit — > 0. 
Such discrepancy is well understood, as in this regime 
the analytical treatment is invalid (as it is based on the 
approximation, that the dynamics of tunneling exceeds 
the rate of electromechanical coupling). 

We have made further studies, considering the situa- 
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tion when the value of Tl is such that the corrections 
scaling with 1/Tl are not negligible, for instance, choos- 
ing Tl = 10w. In this case, we observe a significant 
decrease of the cooling efficiency. In particular, at lower 
values of Q the scheme appears slightly more efficient 
when the resonator is coupled to the left dot rather than 
to the right dot. 

The dynamics is significantly modified, if one chooses 
Tl small and r# comparatively large, see Fig. [5] In 
this regime there is essentially no excess electron in the 
DQD at steady state. While no can be very small, the 
cooling rate is generally very slow, so that one obtains 
relatively large mean phononic numbers at finite values 
of Q. In this slightly better cooling efficiency is 

found when the resonator is coupled to the left dot, see 
Fig.[5](c) and (d). We also note an unexpected behaviour: 
When the resonator is coupled to the right dot, cooling 
is found for negative values of A, as visible in Figs. [5] 
(a) and (b) (hence, when the left dot is at higher energy 
than the right dot). This could be explained in terms 
of interference effects in the electromechanical coupling 
inside the DQD, analogously to the dynamics observed 
in RrfPEI 



C. Discussion 

The scheme here proposed shares several analogies 
with laser cooling schemes of trapped ions, where the 
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FIG. 4: (a) Average excitation number fist and (b) cooling 
rate 7 to t as a function of Tr for 7 P = 0. The solid (dash- 
dotted) curves are found for -y p = (pf p = 10" 5 o;) by maxi- 
mizing the values of 70 as a function of A and T c in Eqs. ( 36 1 
and (35 1. The other parameters are a = O.Icj, Fl = IOOcj, 



= and rip = 50. The dots and the crosses correspond to 
7 P = and 7 P = 10~ j uj, respectively, and are extracted from 
numerical simulation, where the evolution of the resonator 
mean phonon number is calculated by numerical integration 
of master equation ( 10 1, for the same parameters as in the an- 



alytical case. The cooling rate is then determined by fitting 
the curve of the numerical integration with an exponential 
decay, minimizing the sum of the squares of the offsets of the 
numerical points from the exponential fit. The value of fist is 
extracted from the behaviour of the curve at sufficiently long 
evolution times. The dashed line in (b), visible at small Fr, 
indicates the regime in which the analytical results are not 
valid. 



mechanical effects of photon-atom interaction is used in 
order to prepare the motion of trappe d par ticles in the 
ground state of the confining potentia l 45 * 46 !. It is inter- 
esting to draw more explicit analogies, in order to evi- 
dentiate also the differences which arise when electrome- 
chanical effects are used in place of optomechanics. 

A clear analogy between the two systems is found when 
considering the mechanical forces exerted on the har- 
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resonator is coupled to the right dot, whereas in (c),(d) the 
resonator is coupled to the left dot. 



monic oscillator. Indeed, the interaction Hamiltonian in 
Eq. ([4]) can be mapped to the Hamiltonian, describing 
the mechanical effects by photon scattering in the Lamb- 
Dicke regime, by means of the unitary transformatiorP^ 



(41) 



with p — i (<v — a) such that p — pS, and the system 
operators transform as 



5^0)^15 = \0){R\e^ 
S^aS 



a - -\R)(R\. 

UJ 



(42) 



In this reference frame the interaction Hamiltonian is 
given by 

H Te ^ ph = -T c (\L) (R\c-^p + \R) (L\j&) (43) 

which has the same form as the interaction between a 
mode of the electromagnetic field and an atomic dipo- 
lar transition. In particular, in the limit of large Tl 
the state |0), with no excess electrons in the DQD, can 
be adiabatically eliminated, and the DQD dynamics can 
be effectively described in the Hilbert space spanned by 
the state \L) and \R). In this regime, the system is 
analog to a trapped two-level atom with linewidth T' d , 
cooled by a laser with Rabi fr equen cy T c and possessing 
a Lamb-Dicke parameter a/cJ 45 l 46 ( The dephasing rate 



hence plays a similar role of the linewidth of the cooling 
transition, setting the fundamental bound to the final 
temperature one can achieve by means of electromechan- 
ical forces in this setup. Drawing on this analogy, ground 
state cooling, as discussed in this paper, is a form of side- 
band cooling for trapped ions, and it is possible provided 
that the dephasing rate T d < lj. In that case, in this pa- 
per we show how the parameters can be tuned in order 
to achieve the largest ground state occupation. Summa- 
rizing, efficient cooling is obtained for large values of Tl 
and small values of T^j, such that < uj. Large ground 
state occupations are also found when Tr ~ uj. When 
Td > uj, it is generally not possible to reach large ground- 
state occupations, even though cooling can be performed 
which has analogous efficiency as Doppler cooling. 

Differing from photon scattering, in the case of elec- 
tron tunneling we do not have diffusion processes, which 
are otherwise encountered when energy is dissipated by 
spontaneous emission in free space. For this reason, 



Eq. ( 38 ) has the same form found in sideband cooling 
of trapped ions under specific condition s, wh ere diffusion 
due to photon scattering is suppresse d 43 * 47 ! . The same 
equation (for other physical parameters) was derived in 
theoretical treatments of cooling of nanomechanical res- 
onators via photons^HH] where the cooling dynamics was 
mapped to sideband cooling of trapped ions. In our case, 
however, the parameters for which fz-o is minimum do 
not coincide with the ones, at which 70 is maximum, 
and optimal cooling is found as a compromise between 
maximizing the electromechanical cooling rate, 70, and 
minimizing the lower bound to the mean phonon value, 
no, so to effectively counteract the thermalization rate 
due to the coupling with the external environment. 



IV. CONCLUSIONS 

We have presented an extensive analysis of the effi- 
ciency of ground state cooling of the phononic mode of a 
resonator, constituted by a suspended carbon-nanotube 
in a double quantum dot configuration. The analysis 
takes into account dephasing in the current through the 
double dot and the finite Q of the resonator. Ground 
state cooling is possible provided that the dephasing 
rate is smaller than the resonator frequency. Moreover, 
largest efficiency arc obtained when the permanence time 
scale of the excess electron inside the double dot exceeds 
the time scale in which the double dot has no excess 
electron. When the dephasing rate is larger than the res- 
onator frequency, the resonator can be still cooled pro- 
vided that the tunable system parameters are accordingly 
chosen. 

Most parameters of the proposed setup can be tuned. 
The frequency uj of the resonator, for instance, depends 
on the length of the nanotube section that is suspended. 
The frequency difference A and the tunneling rates Tc, 
r^,, and Tji, can be controlled by the external gates^. 
Electron dephasing rate, however, remains to be mea- 
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sured in nanotubes. We remark that the proposed cool- 
ing scheme can be applied to other device layouts, such as 
mechanical resonators electrostatically coupled to fixed 
DQDsP. 
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